Non-equilibrium temperatures in steady-state systems with conserved energy 
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We study a class of non-equilibrium lattice models describing local redistributions of a globally 
conserved quantity, which is interpreted as an energy. A particular subclass can be solved exactly, 
allowing to define a statistical temperature Tth along the same lines as in the equilibrium microcanon- 
ical ensemble. We compute the response function and find that when the fluctuation-dissipation 
relation is linear, the slope of this relation differs from the inverse temperature T^^^ . We argue 
that Tth is physically more relevant than Tfd, since in the steady-state regime, it takes equal values 
in two subsystems of a large isolated system. Finally, a numerical renormalization group procedure 
suggests that all models within the class behave similarly at a coarse-grained level, leading to a new 
parameter which describes the deviation from equilibrium. Quantitative predictions concerning this 
parameter are obtained within a mean-field framework. 

PACS numbers: 05.20.-y, 05.70.Ln, 05.10.Cc 
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I. INTRODUCTION 

The existence and the precise definition of intensive 
thermodynamical parameters in out-of-equilibrium sys- 
tems still remains an open issue. Indeed, the goal of a 
statistical approach for non-equilibrium systems, which 
remains to be constructed, would be to give a well-defined 
meaning to such thermodynamical parameters, and to 
predict their relation with extensive macroscopic vari- 
ables like energy or volume. Accordingly, many attempts 
have been made to define out-of-equilibrium tempera- 
tures in the last decades . 

In the context of glasses, which are non-stationary sys- 
tems with very large relaxation times, effective temper- 
atures have been first introduced as phenomenological 
parameters allowing to account for experimental data 
0, 0, ■ More recently, the notion of effective tempera- 
ture has been given a more fundamental status, being de- 
fined as the inverse of the slope of fluctuation-dissipation 
relations (FDR) in the aging regime This definition 
was guided by the dynamical results obtained within a 
family of mean-field spin glass models '^l. Interestingly, 
such a definition of the effective temperature has been 
shown to satisfy the basic properties expected for a tem- 
5ince then, a lot of numerical simulations 
I El m and experiments CI El El 113, El 
have been conducted to test the validity of this definition 
of temperature in the aging regime of glassy materials. 
Yet, this definition seems not to be always applicable as 
the measured FDR can be non linear. 

Other classes of systems are far from equilibrium not 
due to their slow relaxation towards the equilibrium 
state, but rather because they are subjected to external 
constraints producing fluxes (of particles or energy for 
example) traversing the system, leading to energy dis- 
sipation. As a result, they never reach an equilibrium 
state. Among these systems, one can think of granu- 
lar gases, sheared fluid or certain kinetic spin models, to 
quote only a few of them. 
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Although the usual formalism of equilibrium statisti- 
cal physics does not apply to these systems, it is interest- 
ing to note that the latter sometimes share with equilib- 
rium systems some quantitative properties, like critical 
behavior 0, . To describe the statistical properties 
of such non-equilibrium systems, effective temperatures 
have been deflned either from FDR l2ll l2g | or from 
maximum entropy conditions 24], as originally pro- 
posed by Jaynes j^^. Still, the vahdity of these proce- 
dures remains to be clarified in the context of non glassy 
out-of-equilibrium systems. 

When described in a probabilistic language, a com- 
mon feature of these systems is that they do not obey 
the detailed balance property, considered as a signature 
of equilibrium dynamics. Since the breaking of detailed 
balance plays an important role in non-equilibrium sys- 
tems, it may be useful to distinguish between different 
forms of detailed balance which should not be confused. 
In the literature, the term 'detailed balance' often refers 
to a canonical form which reads: 



W{P\a)e 



-Ea,/T _ 



W{a\P) 



(1) 



where T4/'(/3|a) is the transition rate from state a to state 
/5. This criterion on transition rates ensures that the sta- 
tistical equilibrium reached at large times by the system 
is indeed the canonical equilibrium at temperature T . 

Still, the above approach requires to know the equilib- 
rium distribution before defining the stochastic model. 
On the contrary, one could try to find a stochastic model 
which describes in the best possible way a given com- 
plex hamiltonian system, without knowing a priori the 
equilibrium distribution. Such a stochastic model should 
at least preserve the symmetries of the original hamil- 
tonian system, which are the energy conservation and 
the time-reversal symmetry t ^ —t (additional symme- 
tries -translation, rotation, etc.- must be taken also into 
account when present). Energy conservation is easily im- 
plemented in the stochastic rules by allowing only tran- 
sitions between states with the same energy. On the 
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other side, the time-reversal symmetry in the hamilto- 
nian system can be interpreted in a stochastic language 
as the equality between two opposit transition rates: 
VF(/3|a) = iy(Q;|/?), a property called microcanonical de- 
tailed balance or microreversibility. 

In the context of non-equilibrium systems, one expects 
that the time-reversal symmetry is broken due to the 
presence of fluxes or dissipation. Hence, a simple way to 
define a non-equilibrium system is to consider more gen- 
eral microcanonical forms of detailed balance relations 
such as: 



W{P\a)f^ = W{a\(3)ff, 



(2) 



where is the statistical weight of state a, and with 

Ea = Efj. 

In this paper, we study a class C of non-equilibrium lat- 
tice models describing local redistributions of a globally 
conserved quantity, which is interpreted as an energy. 
A particular subclass Cs satisfies a microcanonical de- 
tailed balance relation of the form (01 , which differs from 
microreversibility and can be solved exactly, allowing to 
define a statistical temperature Tt^ along the same lines 
as in the equilibrium microcanonical ensemble. The re- 
sponse function is computed explicitely, and the FDR is 
found to be linear or non linear depending on the model 
considered -the response can even be non linear with the 
perturbing field. Very interestingly, when the FDR is lin- 
ear, its slope differs from the inverse temperature T^^^, 
which questions the relevance of FDR to define a temper- 
ature in non glassy out-of-equilibrium systems. Finally, 
we implement numerically a functional renormalization 
group procedure to argue that all the models within the 
class C behave at the coarse-grained level as a member of 
the subclass Cs- Predictions about the renormalization 
procedure are also made using mean-field arguments, and 
are quantitatively verified. Note that a short version of 
some aspects of this work has appeared in , and that 
a related model, including kinetic constraints, has also 
been introduced in the context of glassy dynamics . 



II. MODELS AND STEADY-STATE 
PROPERTIES 

A. Definition 

The models we consider in this paper are defined as 
follows. On each site i of a d-dimensional lattice, a real 
variable Xi which can take either positive or negative val- 
ues is introduced. Dynamical rules are defined such that 
the quantity 



N 



E 



51 5(2;*) 



(3) 



is conserved. The function g{x), assumed to be positive 
with continuous derivative, decreases for x < xq and in- 
creases for X > Xq, where xq is an arbitrary given value. 



Without loss of generality, we assume g{xQ) = 0. Quite 
importantly, the steady-state distribution can be com- 
puted with these hypotheses only. However, to clarify 
the presentation, we assume in this section that xq = 
and g{x) is an even function of x. 

It is also necessary to introduce the reciprocal function 
g~^{y), given as the positive root of the equation g{x) = 
y. The dynamics is defined as follows: at each time step, 
a link (j, k) is randomly chosen on the lattice and the 
corresponding variables Xj and Xk are updated so as to 
conserve the energy g{xj) +g{xk) of the link. To be more 
specific, the new values x'j and x'^. are given by 



= ±g-\{l ~ q)S,k) (4) 



with Sjk = g{xj) + g{xk), and g is a random variable 
drawn from a distribution ipil)! assumed to be symmetric 
with respect to q = ^ {0 < q < 1). The new values x'j and 
x'j, are either positive or negative with equal probability, 
and without correlation between the signs. Thus a model 
belonging to this class is characterized by two functions 
g{x) and ip{q)- 



B. Master equation 

The system is described by the distribution P{{xi}, t), 
which gives the probability to be in a configuration {xi} 
at time t. Its evolution is given by the master equation: 

BP r 

-g^{{x.},t) = Jl[dx',W{{x,}\{^'SPiWi},t) (5) 
l[dx',Wi{x'Mx.})Pi{x.},t) 



where is the transition rate from configura- 

tion {xi} to configuration {x[}. The transition rate can 
be decomposed into a sum over the links of the lattice: 

wi{x'M^.})^Y.^MWM^^}) (6) 

U-k) 

where VFjfe({a;i}|{a;i}) accounts for the redistribution 
over a given link {j, k): 



W,k{{x'Mx,}) 



dq^j{q) X (7) 



7 H^j-'^]9 ^{iSjk)] S[x'k-<Jkg ^{{l-q)Sjk)] 



where the variables cTj, ak = il account for the random 
signs appearing in Eq. Q with probability ^ -hence the 
factor J in the above equation. After some algebra, the 
transition rate W({a;-}|{a:i}) can be rewritten as: 



n ^(^' 



(8) 
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where g'{x) denotes the derivative of 5(2;). 

C. Detailed balance and steady-state distribution 

A case of particular interest is the subclass of models 
for which the distribution ip{q) is given by a symmetric 
beta law: 



(9) 



with rj > 0. In this case, the function 'ip{g{x'j)/ Sjk) ap- 
pearing in the transition rates factorizes, if one takes into 
account the delta function. So the transition rate reads 



4r(?7)2 



5{x[ - Xi) 



X (10) 



|g'K-)g'K)l 



g{x'^yg{x'kT-H +5(4) - S.u] 



From this last expression, it can be checked that a de- 
tailed balance relation is satisfied: 



N 



W{.{x[}\{x,})\{[\g'{x,)\g{x,T-'] = (11) 

i=l 

N 

W{{x.}\{x[})X{[\g'{x[)\g{x',y-'] 



i=l 



As a result, the steady-state distribution Pst{{xi}\E), for 
a given value E of the energy, is readily obtained as: 



Pstiix^m = .^l-f[[\g'ix,)\g{x.r-']x (12) 



s(j2gix,)-E^ 



where Z]^{E) is a normalization factor that may be called 
an effective (microcanonical) partition function: 



N 



(13) 



ZNiE) = f ll[dx,\g'{x,)\g{x,y>-'] x 
i=i 

6(j2gix.)-E^ 



An important remark has to be made at this stage: 
Eqs. H12|l and 1)1311 remain formally valid if one slightly 
changes the definition of the model. This can be done 
in two different ways. First, one could consider the case 
where the variables {xi} take only positive values. Then 



one only needs to remove the sum \J2a transi- 
tion rates given in Eq. lO, and Eq. (|12|l is recovered, with 
this time Xi > 0. Second, as mentioned in Sect. Ill Al the 
model can be generalized by assuming that g{x) is not 
an even function. This is particularly useful if one wants 
to include an external field which breaks the +/— sym- 
metry -see Sect. IIII Bl Actually, if g(x) decreases for 
X < xq, and increases for x > Xq, the distribution given 
in Eq. (O also holds jH. 

The function Zn {E) can be computed using a Laplace 
transform. Indeed, it appears rather clearly from 
Eq. (|13ll . by making the change of variable Si — g{xi), 
that Zpj(E) is actually independent of the functional 
form of g{x). One finds 



Zn{E)^knE'^^- 



(14) 



with KN = 2^r(77)^/r(7/7V). The fact that Zn{E) does 
not depend on g{x) is actually not a coincidence, but 
comes from the basic definition of the model given in 
Eq. Q). Indeed, for any function g{x), one could choose 
as the dynamical variables the local energies Si = g{xi), 
and solve the model for e^. Coming back to the variable 
g{xi) at the end of the calculations, the distribution H12(l 
would be recovered. Still, it should not be concluded 
from this that all physical quantities defined in the model 
are independent of g{x). In particular, the response to 
a perturbing field depends strongly on g{x), since the 
field is coupled to x, and not to the energy g(x) -see 
Sect. HlTBl 

An interesting question is also to see under what con- 
ditions microreversibility (to be associated to the equilib- 
rium behavior) can be recovered in this model. Microre- 
versibility holds if \g' {x)\ g{x)^~^ is independent of 
can be seen from Eq. (|12|l . Such a condition can be satis- 
fied only if g(x) is a power law, say g{x) — x^/p, where p 
is an even integer to ensure the regularity of g{x) around 
X ~ 0. The factor 1/p has been added for convenience, 
but is otherwise arbitrary. One then has 



\g'{x)\g{xr-'=p'-^\xr- 



(15) 



Accordingly, microreversibility is recovered for rj = 1/p. 
On the contrary, for i] ^ 1 /p, significant differences with 
the equilibrium behavior are expected. These differences 
may be even stronger if g{x) is not a power law. 



III. NON-EQUILIBRIUM TEMPERATURES 
A. Statistical approach 

1. Microcanonical equilibrium 

In order to define a temperature in this model, one 
can try to follow a procedure similar to that of the mi- 
crocanonical ensemble in equilibrium statistical physics. 
Indeed, one of the main motivations when building the 
present model was to find a model in which a global 
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quantity (the energy) is conserved, so as to 'mimic' in 
some sense a microcanonical situation. Yet, as mentioned 
above, the absence of microreversibihty should yield im- 
portant differences with the latter case. For an equilib- 
rium system in the microcanonical ensemble, tempera- 
ture is introduced in the following way. Considering a 
large system S with fixed energy, one introduces a par- 
tition into two subsystems Si and 52, with energy Eg 
and a number Ng of degrees of freedom {£ = 1,2). These 
two subsystems are no longer isolated, since they can 
mutually exchange energy; the only constraint is that 
El + E2 = E is fixed. The key quantity is then the 
number ^lN^{Ei) of accessible states with energy E^ in 
the subsystem Si; in systems with continuous degrees of 
freedom (like a classical gas for instance), VlNi{Ei) is the 
area of the hypersurface of energy E( in phase space. As- 
suming that both subsystems do not interact except by 
exchanging energy, the number of states of the system S 
compatible with the partition (i?i,i?2) of the energy is 
equal to ntq^{Ei)nif^{E2) ■ But since Ei + E2 is fixed, 
the most probable value E^ is found from the maximum, 
with respect to Ei, of VIjm^{Ei)VIjm^{E — Ei). Taking a 
logarithmic derivative, one finds the usual result: 



dEi 



N2 



dEo 



E-ET 



(16) 



Defining the microcanonical temperature Ti of subsystem 
£ by the relation 



1 



dlnQ 



Ne 



dEo 



(17) 



one sees from Eq. (|16|) that Ti — T2, i.e. that the tem- 
peratures are equal in both subsystems (throughout the 
paper, the Boltzmann constant ks is set to unity). In 
addition, it can also be shown that the common value T 
does not depend on the partition chosen; as a result, T 
is said to characterize the full system S. 



2. 'Microcanonical' stationary state 

Very interestingly, this microcanonical definition of 
temperature can be generalized in a rather straightfor- 
ward way to the present model. Still, it should be noticed 
first that microscopic configurations compatible with the 
given value of the energy are no longer cquiprobablc, as 
seen from the distribution (|12|1 . so that ^Im{E) is no more 
relevant to the problem. But starting again from a par- 
tition into two subsystems as above, one can determine 
the most probable value E^ from the maximum of the 
conditional probability P{Ei\E) that subsystem Si has 
energy Ei given that the total energy is E. Indeed, in 
the equilibrium case, P{Ei\E) reads 

= (18) 

which by derivation with respect to Ei, yields precisely 
the same result as Eq. (|16ll . 



To be more specific, the subsystems are defined in the 
present model as a partition of the lattice, with Ni sites 
in Si and A'2 sites in ^2. The conditional distribution 
P{Ei\E) is then given by: 



P{Ei\E) 



J fldx,Psti{x,}\E)6ij29^^ 

. N 

Y[[dx,\g'ix.,)\g{x.,)'^-' 



1 



Zn{E) 



N 



E 5 




Taking into account the last delta function, the first one 
can be replaced by gi^i) — {E — Ei)), so that 

P{Ei\E) may be written in a compact form as: 



P{Ei\E) 



Zni {El ) {E — El) 
Zn{E) 



(20) 



This result generalizes in a nice way the equilibrium dis- 
tribution Eq. H18(l . since in equilibrium Zjq{E) reduces 
precisely to Q,m{E). The most probable value E^ satis- 
fies 



ainP(£;i|£;) 



dEi 







which yields 



d In Zjsi^ 



dEi 



51nZ 



N2 



dE2 



E~ET 



(21) 



(22) 



So in close analogy with the equilibrium approach, we 
define a temperature T/^ for subsystem Sg through 



1 



dlnZ 



th 



dEf 



(23) 



Then Eq. ^ imphes that T/^ = T^^. 

At this stage, it is important to check that the com- 
mon value Tth of the temperature does not depend on 
the partition chosen. To this aim, we show that Tth can 
be expressed as a function of global quantities character- 
izing the whole system, with no reference to the specific 
partition. 

Let us compute Zpj(E) as a function of Zn-^{Ei) and 
Z]^^{E2). Since dEi P{Ei\E) = 1, one has from 
Eq."|29: 



Zn{E)= [ dEiZMAEi)ZNAE~Ei) 
Jo 



(24) 



We assume the following general scaling form at large Ng 
for ZN,{Ei): 



ZnAEi) = Ai e^p[NeCeisi)] 



(25) 
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with = Ei/Ng (the index ^ = 1,2 labels the subsys- 
tem). This scaling form is demonstrated explicitely in 
Sect. IIII A 31 Using a saddle-point calculation, one ob- 
tains for Zn{E) a relation of the form (with e = E/N): 

Zn{E)^ZnAEI)ZnAE*2)N rd£ie-^''(^)(^i-^^)' 

(26) 

where 6(e) is defined as: 

m = -\[^iC'l{el) + ^2C'i{e--el)] (27) 
with \i = Ni/N. Thus \rZn{E) reads: 



In Zn{E) = In Zn, {E*) + In Zn, {E - E*) - - In b{e) + C 

(28) 

where C does not depend on E. Taking the derivative 
with respect to E yields, using Eq. H23|l : 

ain^_ 1 dE* 1 / dEl\ 1 

In the limit N ^ oo (with s fixed), the last term vanishes, 
whereas Tth has a finite limit due to the scaling form 
Eq. idll), so that 



d\nZ 



N 



dE 



1 

Tth 



(30) 



As a result, Tth can be computed from the global quan- 
tity Zn{E) instead of Zjv^(iJi) or Zm2[E2)^ and is thus 
independent of the partition chosen. This temperature 
characterizes the statistical state of the whole system. 
From Eq. H14|l . the equation of state of the system is: 

E = r]NTth (31) 

In the case of a quadratic energy, i.e. g{x) = ^a;^, it 
has been shown above that the equilibrium behavior is 
recovered for 1] = ^. This result is confirmed by Eq. (|31|l . 
which reduces for = i to the usual form of the energy 
equipartition. On the contrary, for 77 ^ i, a generalized 
form of equipartition holds in the sense that all the sites 
have the same average energy e = E/N (which is not 
surprising given the homogeneity of the system), but this 
average energy per degree of freedom is equal to rjTth 
instead of ^Tth- This point will be discussed in more 
details later on. 

Up to now, we have considered only the 'microcanon- 
ical' (in a generalized sense) distribution Pst{{xi}\E). 
Yet, it would be interesting to introduce also the analo- 
gous of the canonical distribution. To do so, we compute 
the distribution Pcan{{xi}) associated to a small (but still 
macroscopic) subsystem Si of a large isolated system S. 
The degrees of freedom {xi} with i = Ni + 1, . . . , N have 
to be integrated out since they belong to the reservoir. 



One finds for the remaining {xi} (i = 1, . . . , A^^i) the fol- 
lowing distribution: 



Pcaniix.}) = y4^n 9ix^r-' X (32) 



The above integral is nothing but the partition function 
Zn2{E - J2f=i 9{xi)), with N2 = N - Ni, which can be 
expanded to first order as: 

\nZM2 f E -^g(xi) \ = InZM^iE) - — ^9{xi) 

\ i—1 / i—1 

(33) 

assuming that ff(^j) ^ E, which is true as long as 

Ni ^ N. The derivative of In Z^^ [E] has been identified 
with 1/Tth using Eq. (|23l) . up to corrections that vanish 
in the limit Ni/N 0, since E is the total energy rather 
than the energy E2 of the reservoir. Introducing this last 
result into Eq. (|32() . one finally finds 

Pcan{{x^}) = l[\g'{x,)\g{xi}^-^X (34) 



exp 




where = Zn^ (E) /Zn{E) -note that E is the energy 
of the global system which includes the reservoir. This 
'canonical' distribution appears to be useful in order to 
compute the FDR, as discussed below in Sect. IIII Bl 



3. Entropy and thermodynamics 

From Eq. H3U|) , it is tempting to generalize the notion of 
microcanonical entropy through S{E) = In Zf^(E). In- 
deed, this definition is not only an analogy, but as we 
shall see, it can be associated with a time-dependent en- 
tropy which is maximized by the dynamics. To define the 
entropy, one needs first to introduce the probability mea- 
sure PE{{xi},t) restricted to the hypersurface of energy 
E: 

Pi{x,}, t) = Pe{{x,}, t) 5 g{xi) - E^ (35) 
Then the dynamical entropy is defined as: 



N 

SE{t) = - J Y[dx,P({x^},t) In 



PE{{x^},t) 
f{{x^}) 



(36) 



where /({xj) = Hili W ixi)\ 9{xi)'^ ^ Using the mas- 
ter equation ||SJ|, it can be shown that SE{t) is a non- 
decreasing function of time -see Appendix ^ As a re- 
sult, SEit) is maximal in the stationary state, and the 
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corresponding value S{E) is given by: 



SiE) 



N ^ 

j X{dx,P{{x.^\E) In — 



\nZN{E) 



(37) 



which matches exactly the definition proposed above on 
the basis of Eq. (ppi . 

Using Eq. H14|l . one can compute S{E) and check ex- 
plicitely that the entropy per site S{E)/N becomes in 
the thermodynamic limit a well-defined function of 
the energy density e = E/N. The entropy S{E) reads: 



S{E) ^ N 



In 2r(r/) - — In r(?7iV) + In S 



As lnr(a::) w xlnx — x for large x, one finds: 
In r(r/iV) w 77(ln 77 - 1) + ry In TV 

which allows to write S{E) — NC{e) with: 

C(e) = 7?ln£ + ln2r(77) - 7?(lnr/ - 1) 



(38) 



(39) 



(40) 



On the other hand, the equilibrium thermodynamic 
formalism is most often formulated in terms of the canon- 
ical ensemble. In the present model, since a canonical 
distribution has been derived, it may also be possible to 
define an equivalent of the canonical thermodynamic for- 
malism. Indeed, from Eq. (|34|l . one can easily see that 
the average energy (E) is given by 



ainZ^v"" 
dp 



(41) 



where /3 = T^^^ is the inverse temperature. A generalized 
free energy F{Tth) is also naturally introduced through 



can 
N 



The generalized partition function Z 
computed, as it is factorized: 



can 
N 



(42) 

can be easily 



yean 



dx \g'{x)\g{x) 



N 



which leads to 

= (2r(r;)T,'',)^ 
So the free-energy is given by 

F^-NTth [\n2T{ri)+f^\nTth] 



(43) 



(44) 



(45) 



In equilibrium, the entropy 5* is related to the free energy 
F through 



OF 
dT 



-S 



(46) 



This relation is also satisfied within the present model: 
OF 



dTt 



th 



N[r]{\nTth + l)+lr^2T{ri)] 
NCie) (47) 



where the last equality is obtained by using the equation 
of state Tth = 'S/ij, and comparing with Eq. (|40|l . 



B. Fluctuation-dissipation relations 

As recalled in the introduction, temperatures are usu- 
ally defined in out-of-equilibrium systems as the inverse 
slope of the FDR, when this relation is linear. This ap- 
proach has been shown to be physically meaningful in 
the context of glassy models in the aging regime 6] . In 
this case, the long time slope of the FDR gives an effec- 
tive temperature which differs from the heat bath tem- 
perature. Still, for non-equilibrium steady-state systems 
which are not glassy, no justification has been proposed 
to show that the inverse slope of the FDR satisfies the ba- 
sic properties expected for a temperature. For instance, 
one expects a temperature to take equal values in two 
subsystems of a large system, when the stationary state 
has been reached. The present model thus allows to test 
explicitely the validity of the FDR definition of temper- 
ature. 

A natural observable to consider in this model is 



N 



M(t) = ^x.(i) 



(48) 



The steady-state correlation function C{t) of the system 
is then defined as the normalized autocorrelation of the 
observable M (t) between time t — and t: 

at) = l((M(t) - {M}) (Af(0) - (M))) (49) 

where the brackets (...) denote an average over all pos- 
sible trajectories of the system. Calculations are easier 
using the canonical distribution Pcan{{xi}); since this 
distribution is factorized, the random variables Xi and 
Xj are independent if z 7^ j, so that C{t) reduces to 



C{t)^{{x{t)-{x)){x{0)^{x))) 



(50) 



where x stands for any of the variables Xi -all sites have 
the same average values. 

The aim of the FDR is to relate correlation and re- 
sponse of a given observable. One thus needs to introduce 
a perturbation which generates variations of x so that a 
response could be defined. A simple way to perturb the 
system is to add to the energy a linear term proportional 
to an external field h: one then replaces E by E^ defined 
as: 



N N 

Eh ^^9hix^) = 

1=1 i=l 



g{xi) - hx., + Ch (51) 
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Without loss of generality, the new function gh{x) is 
shifted by a constant Ch so that the minimum value of 
gh{x) remains equal to 0. If the second derivative g"{0) 
does not vanish, Ch is given to leading order in h by: 



Ch 



2g"(0) 



(52) 



In order to define the response function, one assumes that 
the system is subjected to a field h ^ for t < 0, and that 
it has reached a steady state. Then at time t = 0, the 
field h is switched off. The (time-dependent) response is 
defined for t > Q through: 



x{t) 



d_ 

mi 



h=o \ N 




(53) 



where the index h on the brackets indicate that the av- 
erage is taken over the dynamics in presence of the field 
h. The observable {N~-^ J2i^i{^))h can be computed as: 




TV 



l[dx.dx',G'i{{x.}\{x[}) X (54) 



where {{xi}\{x'^}) is the zero-field Green function, 
i.e. the probability for the system to be in a configu- 
ration {xi} at time t, given that it was in a configuration 
{x[} at time t = 0, in the absence of field. The response 
function x(t) is obtained by taking the derivative of the 
above equation with respect to h, at h = 0: 



N 

it) = J l[dx,dx',GU{Xr}\{x',}) 



(55) 



.(K},0)^^^(K},0) 



N ^ 



i=i 



The canonical distribution Pcani{xi} , h) in the presence 
of field takes the same form as Eq. H34|l . simply replac- 
ing g{x) by gh{x). Thus one finds for the logarithmic 
derivative of Pcan{{xi}, h): 



dh 



K{^a,o) 



N 

E 



dh 
( 



/i=0 

I 



(56) 



^ \Tth g'{xi) 



9{xi) 



Note that dch/dh = at /i = 0, due to the regularity of 
g{x). The derivative of the partition function yields: 



dh 



h=0 



I N 

E' 



(57) 



where cj, stands for 



9'ix^) 



+ (^-1) 



9ix^) 



(58) 



Replacing the expression (|56|l in Eq. H55|l . one finally 
finds, using the factorization of the canonical distribu- 
tion: 



x{t) 



1 

Tth 



C{t) - {{x{t) - (x))(^(0) - (^))) (59) 



where indices are omitted just as in Eq. (|50|l . Com- 
pared to the usual form of FDR, an additional term ap- 
pears which corresponds to the correlation of the vari- 
ables X and Lo. In general, this new correlation function 
is not proportional to C{t), so that a parametric plot 
of xit) versus C{t), usually referred to as a fluctuation- 
dissipation plot, would be non linear. 

Yet, in the case where g{x) is an even function of x, 
some important simplifications occur. On the one hand, 
the average values of x and to vanish. On the other 
hand, the correlation x{t)uj{0) becomes proportional to 
the 'hopping correlation function' ^{t), defined as 



/ 1 



(60) 



The variables 4'iit) ^-^e history dependent random vari- 
ables, which are equal to 1 if no redistribution involv- 
ing site i occured between t = and t, and are equal 
to otherwise. The proportionality of both correlation 
functions can be understood as follows: if there was a 
redistribution on site i between and t, Xi{t) becomes 
fully decorrelated from Ui, due to the fact that the sign 
of Xi{t) is chosen at random, and that the average values 
(x) and (lo) vanish for an even g{x). On the contrary, if 
no redistribution occured, Xi{t)u}i{0) = Xi{0)uji{0). The 
same reasoning also holds for C{t), so that one has: 

C(t) = (x^) $(t), {x{t)uj{Q)) = (xuj) $(t) (61) 

As a result, the FDR can be expressed, in the case of an 
even function g{x), as 



Cit) 



(62) 



So the FDR is indeed linear in this case, and one can de- 
fine an effective temperature Tpu from the inverse slope 
of this relation. This yields: 



1 



FD 



1 

Tth 



(XLO) 



(63) 



Still, as long as {xuj) ^ 0, the temperature Tfd differs 
from the temperature Tth defined above from statistical 
considerations -a more detailed discussion on this point 
is given below in Sect. IIIlUl 

Even though the two temperatures are not equal, one 
can wonder whether they are proportional, in the sense 
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that the ratio Tth/Tpo would be independent of Tth- 
From Eq. H63|l . one has: 

^ = 1-^^^^^ (64) 

where we have used the relations Tth ~ ^/v and e ~ 
{g{x)). The correlation {xu) can be written in a more 
explicit form as 

(xc.) = /^ + (r7-l)^\ (65) 
\g'{x) g{x)/ 

From Eqs. (|64|l and (|65() . it appears that the ratio 
Tth/TpD generally depends on Tth, since the average (. . .) 
is done with the one-site distribution which is a function 
of temperature. 

Now in the particular case where g{x) is a power law, 
namely g{x) — x^/p (with p an even integer), Eq. H64fl 
actually simplifies to 

TFD--[2+piv-l)]Tth (66) 

Note that for 2+p{r] — 1) < 0, the above equation would 
lead to a negative TpD, i.e. a negative response xit) to 
the perturbation h, which is rather counterintuitive. Ac- 
tually, x{t) does not become negative in this case but 
diverges and Eq. H66|l is no longer valid, indicating the 
breakdown of linear response -the response is then non 
linear with h even for h 0. This may be seen from the 
correlation {xuj), which can be written: 

/>oo 

{xuj)^A dxa;i+P("-i)e-^''/f^*'' (67) 
Jo 

where the constant A depends onp and ry. If l+p^rj—l) < 
— 1 (i.e. the same condition as above), the integral di- 
verges at its lower bound, and xit) becomes infinite. To 
keep the susceptibility finite, one needs to consider val- 
ues of 7] such that rj > 1 — 2/p. It is interesting to note 
that as soon as p > 2, the equilibrium value rj — 1/p 
does not satisfy the above inequality, so that the equilib- 
rium response is non linear in this case. This is somehow 
reminiscent of the Landau theory for phase transitions, 
in which the magnetization (m) becomes non linear with 
the magnetic field at the critical point, where the term 
in in the expansion of the free-energy vanishes. 

Finally, considering the specific case g{x) = as in 
[26|. the above restriction disappears since 1 — 2/p = 
for p = 2. The temperature Tpo is then defined for all 
ry > [s^l- Using Eqs. lpTT|l and one can write Tfd 
in a very simple form which docs not depend on rj: 

Tfd = 2e (68) 

where e is the energy density e — ^{x'^)- 

To sum up, several different cases have to be distin- 
guished. For general regular functions g{x) with g{x) ~ 
x^ for a; — > 0, where p > 2 an even integer, the response 
is non linear with the field ft, if 77 < 1 — 2/p. Other- 
wise, the response is linear and the susceptibility x(i) 



can be defined. In this case, the FDR (or equivalently, 
the fluctuation-dissipation plot) is generically non lin- 
ear. Now, several additional assumptions on g{x) can be 
made: if g{x) is even, the FDR is linear, leading to the 
definition of Tfd as the inverse slope of the FDR; yet, 
TpD is a priori not proportional to Tth- Besides, if g{x) 
is a power law (and if the response is linear), then Tfd 
becomes proportional to Tth- The equality Tfd = Tth 
is recovered only for p = 2 and rj — i, i.e. when linear 
response and microreversibility hold. 



C. Physical relevance of the different temperatures 

In the preceding sections, two different temperatures 
have been introduced: a first one (Tth) from statistical 
considerations, and a second one {Tfd) from a FDR. 
These two temperatures do not only have different defi- 
nitions, but they also take different values, as seen from 
Eq. (|64|l . In this section, we wish to compare the physical 
relevance of these two definitions, and see whether or not 
both of them satisfy the basic properties expected for a 
temperature. 



1. Inhomogeneous version of the model 

Considering a homogeneous system as we have done up 
to now, it is clear that if Tth takes the same value in two 
subsystems, so does Tfd since the two temperatures are 
related through Eq. Indeed, if g{x) = x^/p these 

two temperatures are proportional according to Eq. H66|l . 
so that they may be considered to be identical up to 
a redefinition of the temperature scale. As a result, it 
seems not to be possible to discriminate between these 
two definitions within the present model. 

Actually, this apparent equivalence of both tempera- 
tures comes from the fact that the parameter rj is the 
same throughout the system. So one could try to pro- 
pose a generalization of the model in which rj would not 
be constant, still keeping the model tractable. This can 
be realized in the following way. Introducing on each 
site i a parameter 77^ > 0, we define on each link {j,k) a 
distribution ipjk (q) through 

^M-^rK^,<i'''-'i^-'iy''-' (69) 
r(f?j)r(?7fc) 

The redistribution rules are assumed to keep the same 
form as in Eq. Q. Yet, links (j, k) now need to be ori- 
ented since ipjkil) is no longer symmetric, so that the 
fraction q is attributed to site j, whereas 1 — g is at- 
tributed to site k, precisely as in Eq. Q). 

Note however that even though the redistribution pro- 
cess is locally biased if rjj ^ rjk , there is no global energy 
flux in the system since the form (|59|l has been chosen 
to preserve the detailed balance relation. As a result, 
the steady-state distribution can be computed exactly 
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for any set of variables {jyi}. To siraplify the discussion, 
we restrict the results presented here to the simple case 
g{x) — ^ but generalization to other functions g{x) 
are rather straightforward. In this case, the 'microcanon- 
ical' distribution P{{xi]\E) takes essentially the same 
form as previously: 

(70) 

Following the same reasoning as above, one can define 
both and Tpjj in this generalized model. In particu- 
lar, the temperature T^^^ is defined from the conditional 
probability P{Ei\E) as in Eq. (|23|l . Considering again a 
partition of a large isolated system into two subsystems 
Si and 52, one finds for the subsystem Si 

tL = T'fd = (71) 

where (77) £ is the average of rji over the subsystem S^: 

If one chooses the set of variables {rji} such that (7/)i ^ 
(77)2, the equality T^/^ — Tfj^, which is true from the very 
definition of T/^ -see Eq. (I22|) ~ implies £1 7^ £2- Conse- 
quently, equipartition of energy breaks down, and from 
Eq. (|7TT) one has T^^ ^ T'fd ■ fluctuation-dissipation 
temperature does not take equal values in two subsystems 

m- 

This last point is indeed reminiscent of recent numer- 
ical results reported in the context of binary granular 
gases |22:], where the temperature associated to each 
species of grains from a FDR does not equilibrate. These 
results indicate that for non glassy systems, the temper- 
ature defined from FDR does not fulfill the basic prop- 
erties required for a temperature, as the equality of the 
temperatures of subsystems when a steady state has been 
reached. On the contrary, the temperature Tth defined 
from statistical considerations satisfies this property, and 
may thus be given a more fundamental status. 

Finally, it should be noticed that the relation Tth —^/v 
indicates that the temperature Tth is not simply a mea- 
sure of the average energy, but also takes into account 
the fluctuations of energy. Indeed, a large value of 77 
corresponds on the one hand to a low value of the tem- 
perature, and on the other hand to a sharp distribution 
ipiQ), which in turn leads to small energy fiuctuations in 
the system, as can be seen for instance from the canonical 
distribution given in Eq. H34|l . 

2. How to define a thermometer? 

Once a temperature has been formally defined in a sys- 
tem, a very important issue is to be able to measure it, at 



least within a conceptual experiment. This question is in 
general highly non trivial for out-of-equilibrium systems. 
In the context of glassy systems for instance, it has been 
proposed to use a simple harmonic oscillator connected 
to the system as a thermometer Still, in order to 
measure a temperature associated to a given time scale 
T (assumed to be large with respect to the microscopic 
time scale tq), one must use an harmonic oscillator with 
a characteristic time scale of the order of r. In this case, 
the temperature is obtained through the usual relation 
Eosc = ^T, where Eosc is the average kinetic energy of 
the oscillator. For glassy systems, this temperature has 
also been shown to identify with the temperature defined 
from FDR 0. Besides, a numerical realization of such 
a thermometer has been proposed by using a brownian 
particle with a mass much larger than the other particles, 
in a glassy Lennard- Jones mixture under shear 28] . Such 
a definition of temperature is also consistent with the so- 
called 'granular temperature', defined as 2/d times the 
average kinetic energy of the grains 29] {d is the space 
dimension). 

Interestingly, in the present model which is not glassy, 
a somewhat analogous procedure would be to connect a 
new site to the system, and make it interact with the 
other sites using the current kinetic rules of the model; 
this new site would play the role of a thermometer. As- 
suming again g{x) = ■^x'^, the temperature read off from 
the average energy of the thermometer is precisely Tp d ■ 
At first sight, this seems to be in contradiction with 
the above discussion in which we argued that Tth was 
the physically relevant temperature. The paradox comes 
from the fact that we used without justifying it the rela- 
tion Eosc = ■^T to define the temperature T of the ther- 
mometer as a function of the measurable quantity Sosc- 
Accordingly, such a definition does not ensure that T is 
the temperature of the system. 

One of the most important properties of Tth is precisely 
that it takes equal values within subsystems in contact. 
Actually, to obtain Tth, one needs to know the equation of 
state of the thermometer, that relates measurable quan- 
tities like the average energy Eosc to the temperature Tth- 
Indeed, the fact that it is necessary to know the equation 
of state of the thermometer in order to measure the tem- 
perature is not a specificity of non-equilibrium states, but 
is also true in equilibrium situations, in which one must 
know for instance the relation between the height of a liq- 
uid in a vertical pipe and the temperature of this liquid. 
In the same way, the relation Sosc = jT invoked above is 
not obvious in itself, but results from equilibrium statis- 
tical mechanics. As a result, there is no clear reason why 
this last relation should hold for generic non-equilibrium 
situations. 

Yet, an important point must be mentioned at this 
stage. One of the specificity of non-equilibrium states is 
that there is not a unique way to define a thermal contact 
between two systems. In equilibrium, it is usually enough 
to consider the weak interaction limit in which the energy 
associated with the interaction process is very small com- 
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pared to the other energies involved. On the contrary, 
for non-equiUbrium systems, the conservation of energy 
is not sufficient, since the dynamics can be much richer, 
as illustrated by the presence of the parameter ij in the 
present model. Hereabove, we assumed that the new site 
used as a thermometer was driven by the same dynam- 
ical rules as the system it is in contact with. Yet, in 
practical situation, one would rather use a thermometer 
with a known equation of state to measure the tempera- 
ture of another system for which the equation of state is 
unknown. As a consequence, the dynamics of the ther- 
mometer is expected in general to be different from that 
of the system. Determining the properties that a ther- 
mometer has to satisfy in order to measure correctly the 
temperature thus remains an open question. 



IV. RENORMALIZATION APPROACH 



A. Breaking of detailed balance 



1.5 



0.5 



^ I • I • I 

^' / 'A 

'A 
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«^ square box 
o o beta, T] = 1 
+ + beta, T] = 2 
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FIG. 1: Ratio R{ea) of forward and backward probabilities of 
a transition path; R{£a) 7^ 1 indicates a breaking of detailed 
balance. Distributions ipil) used are the sine-like distribution 
(a) and the 'square box' one (o) -see text for details. Results 
for beta tpiq) with rj — 1 (o) and r] = 2 f-l-j are also presented 
for comparison, showing as expected that detailed balance is 
satished for these distributions. 



If V'('z) is different from a beta law, no simple detailed 
balance relation has been found in this model. In the ab- 
sence of such a relation, it is rather hopeless to find the 
stationary distribution Pst{{xi}\E), even though some 
sophisticated algebraic methods have proven to be effi- 
cient in some cases "s^, "si]. Yet, the fact that we were 
not able to find a detailed balance relation in the model is 
not a proof that the relation does not exist. As a result, 
it appears useful to test numerically the existence of non 
zero probability fluxes even in steady state, which would 
clearly demonstrate the absence of detailed balance. 

As discussed in Sect. Ill (Jl the steady-state distribution 
can be fully determined in terms of the dynamics of the 
local energy = g{xi). In the following, we thus use 
these variables Si as the dynamical variables. The dy- 
namics of Si is the same as that of the variables Xi if one 
considers the case g{xi) — Xi, restricting Xi to be positive. 
The detailed balance property is checked by measuring 
with numerical simulations the probability Pati^a, £b, <5e) 
to observe on a given site i a direct transition from a 
value Ei G [ca, £a + 5e\ to a new value e[ G [sf,, Sb + Se] , as 
well as the reverse probability pbai^b, £a, Ss) to go from 
the interval [eb,£b + Se] to the interval [ea,Sa + 6e]. These 
probabilities are actually obtained by averaging over all 
sites i. One then computes the ratio 

n Pab{£a,£b,Se) 

K = -. ^ (Yd) 

Pba(£b, £a,oe) 

which becomes independent of Se in the limit of small Se. 
Besides, a simple parametrization is to set Sb = £a + A, 
and to compute i? as a function of Ea for a fixed value 
of A. Fig. n] presents the numerical results obtained for 
R{ea) with A = E/N, using distributions 'ip{q) which 
differ significantly from beta laws as the sine-like distri- 
bution ^p{q) = ^1 sin(27rg)|, and the 'square box' one, 
^p{q) = 2 tor J < q < J, and ■i/'('Z) = otherwise. 



Beta laws are also shown for comparison. As expected, 
R{ea) = 1 for beta laws, whereas R{£a) 7^ 1 for other 
distributions, showing that detailed balance is broken in 
this case. 



B. Numerical renormalization procedure 

Even though detailed balance is broken microscopically 
when Tp{q) is different from a beta law, one can wonder 
whether the macroscopic properties of the model differ 
significantly or not from that in the presence of detailed 
balance. Indeed, some studies have shown that 

a weak breaking of detailed balance does not infiuence 
the critical properties of particular classes of spin mod- 
els. In the present model, numerical simulations sug- 
gest that even for distributions ip(q) with a behavior 
far from beta laws, no spatial correlations appear within 
two-points functions. Note that this result is also consis- 
tent with the vanishing of two-point correlations in the 
'q- model' for granular matter |32| . which presents some 
formal similarities (although in a different spirit), but 
also important differences, with the present model. In 
particular, the q-model is static, and the role played by 
time here corresponds to the vertical space direction. In 
addition, the dynamics of the q-model is equivalent to 
a synchronous dynamics, and the conserved quantity is 
linear since it represents the vertical component of forces 
between grains. 

In order to test whether macroscopic properties are in- 
fluenced or not by the breaking of detailed balance at the 
microscopic level, one can try to use a renormalization 
group approach. Even though such an approach might 
not seem natural in a context where no diverging length 
scale appears, this is actually a standard way to compute 
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the effective dynamics at a coarse-grained level. Since no 
analytical solution is available for ip{q) different from a 
beta law, one has to resort to numerical simulations. 

To this aim, the following renormalization procedure is 
introduced. The d-dimensional lattice is divided into cells 
(or blocks) of linear size L, and the effective dynamics 
between cells is measured from numerical simulations of 
the microscopic dynamics. To be more specific, when 
running the microscopic dynamics, one has to choose at 
random a link of the lattice at each time step, and to 
redistribute the energy over the link. If both sites of this 
link belong to the same block, then the redistribution is 
only an intra-block dynamics, and corresponds precisely 
to the degrees of freedom that have to be integrated out 
by the renormalization procedure. As a result, nothing 
is recorded during this particular process. 

On the contrary, if the chosen link lies between two 
different cells, then the process is considered as a redis- 
tribution between blocks, and the effective fraction qj^ 
of energy redistributed is computed. Having chosen an 
orientation of the lattice, one can label for instance by 
1 and 2 the two blocks involved in the process. Clearly, 
the total energy of these two blocks is conserved during 
this process. One thus computes the energy and E'^ 
of each block, and defines the effective fraction qu as the 
ratio: 

To obtain the renormalized energy, one should actually 
divide Ei, by the size of the block (so that the energy 
density is conserved), but this is not essential here since 
we consider only energy ratios. The histogram of the 
values of qn obtained when running the microscopic dy- 
namics is recorded, which gives the renormalized distri- 
bution ipL{q)- One would like to test if for large values 
of L, detailed balance is recovered, which would mean 
that the distribution ipL{q) converges (in some sense to 
be specified) towards a beta law. As usual with renormal- 
ization procedures, the correct way to obtain large block 
sizes is not to consider large blocks from the beginning, 
but instead to start from small blocks and to iterate the 
procedure until the desired size is reached. 

As a result, we started from cells of size L = 2 and 
computed successively ip2iq), ip4iq), ''Psiq)^ etc., by ap- 
plying recursively the same procedure with a microscopic 
dynamics defined by the renormalized iphiq) obtained at 
the step before. Numerical results obtained starting from 
an initial distribution -0(g) — ^| sin(27rg)| are shown on 
Fig. m for space dimensions d—1 and d — 2. For L > 4, 
the resulting distributions iphiq) can be very well fitted 
by beta laws, i.e. by a test distribution ^test('z) of the 
form 

V'test(g)-^!|^[<z(l-g)r-' (75) 

with only one free parameter rj^. This parameter ry^ is 
an increasing function of L, which can be easily under- 
stood from the fact that increasing the size of the blocks 




q 



FIG. 2; Renormalized distribution tpL{q) for increasing sizes 
L, in dimension d = 1. Full lines correspond to one-parameter 
fits with beta distributions. Inset: parameter tjl from the fit 
plotted as a function of L'^ for d = 1 (4-) and d = 2 (o); 
dashed line is the mean-Eeld prediction given in Eq. I93fl . 



reduces the fluctuations of the energy from one block 
to another. So if one lets the size L go to infinity, the 
distribution iphiq) eventually converges to a Dirac delta 
function centered on g = i. This means that the beta 
laws found from fitting the data are to be understood 
as pre-asymptotic distributions rather than as true limit 
distributions. 

Very interestingly, the fitting parameter rjh is found to 
be linear with L'^, as seen in the inset of Fig. [3 This 
behavior can be interpreted in the following way, as- 
suming that the initial distribution ipiq) is a beta law 
with parameter 77. As seen from the calculations done in 
Sect. Ill Cl the local distribution of the energy e-i — gixi) 
is given by a gamma law of exponent 77 and scale param- 
eter (3 = l/Tth-. 

= j^er'e-''^- (76) 

The EiS are independent random variables, so that the 
block energies E'b, defined as: 

^ e^ (77) 

block 

are distributed according to beta laws with exponent 
TjL''', where L'^ is the number of sites within a block. 
Then taking the ratio q^, = E^/iE^ + E'^), one obtains 
for qif a beta distribution of parameter rjL'^, as is well- 
known from the properties of gamma laws. 

So starting from a beta law for ipiq), the above analyt- 
ical argument shows that beta laws are again obtained 
from the renormalization procedure, with a parameter 
rjh linear in L'^. Interestingly, the coefficient of propor- 
tionality is precisely the parameter rj of the microscopic 
law ipiq). So when starting from an arbitrary distribu- 
tion ijj{q), it is natural to define an effective parameter 
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r]e from the fitting parameter 77^ as: 



r]L_ 



(78) 



One can then interpret rjf, as the parameter of the micro- 
scopic beta law which would give the same macroscopic 
behavior of the system as the initial distribution 'ip{q). 



C. Mean-field predictions 

In this section, we aim to predict within a mean-field 
framework the effective exponent rje introduced above, 
for an arbitrary distribution ^{q). In a mean-field de- 
scription, one assumes that the two-site steady-state dis- 
tribution P2 (2^1, 2^2) can be factorized as a product of 
one-site distributions: 



P2(xuX2)=Pl{xi)Pi{x2) 



(79) 



This assumption is valid if i}}{q) is a beta law, as can be 
seen from Eq. (PU- For more general ^{q), it remains a 
priori only an approximation. In order to deal with the 
renormalization procedure, it is more convenient to work 
with the distribution p{ei) of the local energy = g{xi), 
rather than with Pi{xi). In terms of the variables {ci}, 
the redistribution rules read: 



q{ej+£k), 



(1-g) {ej+eu) 



(80) 



Numerical simulations show that after a sufficient coarse- 
graining by the renormalization procedure, the renormal- 
ized distribution ipLiq) becomes a beta law with parame- 
ter rji^ — 7]e L'^- The associated renormalized distribution 
of the block energies Eh is then a gamma law with expo- 
nent r]L and scale parameter I3l- 



PLi^b) - ^^^^^ e 



(81) 



The exponent rj^ can be determined from the first and 
second moments of the distribution pL{Eb). Indeed, one 
finds an average value {E^) = t^l/Pl, and a variance 
Var(^6) = r]L/f3l, with Var(£;6) = {E^) - {Ebf . As a 
result, 77L is given by: 



{E,? 



{ED 



{E,Y 



(82) 



If the initial distribution p{£i) is factorized, the block 
energies Ef, are sums of independent random variables - 
see Eq. (|77|l . So the average value and the variance of E^ 
are simply the sums of the average and variance of the 
variables Si. 

(Eb) = (e) i"*, Var(£;6) = Var(e) L"* (83) 

From Eq. (|S^ . the effective exponent rje — Vl/L'^ is thus 
found to be: 



Ve 



Var(e) 



(84) 



So if we know the two first moments of the distribution 
pie), we are able to compute rje. 

To obtain these moments for an arbitrary ipil): we 
use the following steady-state master equation for the 
distribution pie): 



pie) 



deipiei) / de^pie^) x 
JO 
1 

rfgV'('?)'5(£-9(£i +£2)) (85) 





This equation can be considered as describing the redis- 
tribution process over an isolated single link. Yet, it can 
also be derived from a mean-field version of the model, 
in which redistributions can occur over any pair of sites 
of the system -see Appendix ^ 

Introducing the Laplace transform p(s) defined as 



POO 

Pis) = / ii£e~*^p(e) 
Jo 

one can rewrite Eq. H85(l as 

/>! /'OO poo 

Pis) ^ / dqipiq) / deipiei) / de2pie2)e 
Jo Jo Jo 

. 

The integrals over ei and £2 can be factorized into a 
product of Laplace transforms: 



(86) 



;sg(£i+e2) 



Pis) 



dqil^iq)Piqsy' 



(88) 



From the last equation, the successive moments of pie) 
can be obtained, since they are given by the derivatives 
of Pis) in s = 0: 



dp 

- ^Ts 



s=0 



fP 
d^ 



s=0 



(89) 



Note also that by definition, p(0) = 1. Taking the first 
derivative of Eq. H88|) in s = 0, one recovers that (q) = 
i . More interestingly, the second derivative of Eq. H88() 
yields: 



ds^ 



= 2 / dqq^^iq) 
s=0 Jo 



[dp 


)'. 


d^p 




[ds 


s=0 J 


" d^ 


s=0 



(90) 

In terms of moments, the last equation reads: 

(£2) 2{q^) [(£)2 + {eD] (91) 

To compute ije, we only need the ratio (£)^/(£^), which 
is easily found from the preceding equation: 



(£)2 



2{q' 



l-2(q2) 

Taking into account that (q) = i, r/^ is found to be: 



(92) 



Ve = 



1 

8Var(g) 



(93) 
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D. Analytical arguments 

To conclude this section dedicated to renormalization 
group approaches, we wish to give a heuristic analyti- 
cal argument that may help to understand the numeri- 
cal results presented on Fig. |21 As explained above, the 
renormalization can be worked out exactly in the case 
where is a beta law. The numerical procedure shows 
that other distributions ipiq) converge to beta laws un- 
der renormalization. From an analytical point of view, 
it is more convenient to work with the distribution p{e) 
of the local energy, rather than with ip{q). A beta ip{q) 
is associated to a gamma law for p{e) so that it would 
be interesting to check analytically whether an arbitrary 
p{e) converges to a gamma law under renormalization. 
Note that an implicit assumption here is that the iV-site 
energy distribution is factorized, in a mean-field spirit. 

A general calculation for an arbitrary initial distribu- 
tion p{e) is in fact highly non trivial. We thus restrict 
the following calculations to an initial p{£) which differs 
only slightly from a gamma law: 



p{£) ^p^(£) + \5p{e) 



(94) 



where A <IC 1 is an arbitrarily small parameter, and p-i{£) 
is a gamma distribution similar to that used in Eq. H76|l . 
Since the renormalization conserves the average energy, 
p{£) and P'y(e) must have the same average value £ so as 
to become equivalent after renormalization. Taking also 
into account the normalization condition, 5p{e) has to 
satisfy 

/•oc />oo 

de(5p(£) = 0, / d£e5p{£)^Q. (95) 
Jo 

Let M = L'^ be the number of sites in a block. The 
renormalized energy £/j is given by 



(96) 



zGblock 



The distribution of Pi(£_r) is more easily obtained using 
a Laplace transform: 



Pi(s) =p(s/M) 



M 



Obviously, a fixed point for this equation is p{s) — e^"^, 
which leads to p{£) = 6{e — e). The aim of the present 
calculation is to see whether p{e) and p-yis) converge 'in 
the same way' or not toward the delta distribution. 

Replacing Eq. (|94|l into Eq. H97|) and expanding up to 
first order in A, one has: 



SMs)^Mp,{^y'-'6p{^) 



(98) 



Iterating K times the renormalization procedure, one 
gets: 



K-l 



(99) 



The renormalized gcimiricL distribution p-y^n (s) obtained 
after n iterations is given by 



P7,»(s) 



S£ 



-77 A/" 



(100) 



Then Eq. (|99|1 can then be rewritten: 



6pk{s) = M^ (i + -^j 



Using the relation 



lim 1 



se \ 



one ends up with 

Spk{s) « M^e-'^6p 



(101) 
(102) 

(103) 



Expanding 6p{s) in power of s for s ^ 0, one has 

Sp{s) = g2S^ + 0{s^) (104) 

since the terms of order and 1 vanish due to Eq. H95|l . 
This yields: 



SpKis) 



se92S 



(105) 



which goes to when if ^ oo as expected. Yet, this is 
not enough to show that p{e) and Pjis) converge 'in the 
same way' toward the distribution S{e — e). To do so, one 
has to show that 5pk{s) goes to more rapidly than the 
'distance' between p^^k{s) and the infinite K limit 



Pl,oo{s) 



(106) 



A way to quantify this 'distance' is to introduce the quan- 
tity: 

Dk ^ / ds\pj^K{s) - Pj,oo{s)\ (107) 
Jo 

which can be shown easily to take the asymptotic form: 

1 



Dk 



(97) The convergence criterion can be written as 



lim 



Spk{s) 



D 







K 



(108) 



(109) 



This requires that 52 = in the expansion of Spa(s), 
which implies that the distributions p{e) and Pjis) have 
the same variance cr^ = tr^. Such a condition is actually 
natural, as the variance becomes aj^ = a'^/M^ under 
renormalization. If the two distributions take the same 
form after renormalization, they should have in partic- 
ular the same variance cr^ = cr^^, and one recovers 



rt=0 



Obviously, the above arguments are not fully rigorous, 
and remain somehow at a heuristic level, but they al- 
ready give some insights on the mechanisms leading to 
the convergence process observed numerically. 
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CONCLUSION 



supported by the Swiss National Science Foundation. 



The class of models studied in the present paper is a 
very interesting example in which one can define a mean- 
ingful temperature Tth from the conditional energy distri- 
bution of two subsystems, a procedure similar to the one 
used in the equilibrium microcanonical ensemble. These 
models exhibit a rich behavior which includes linear as 
well as non linear response to a perturbation, and linear 
or non linear fluctuation-dissipation relations when the 
response is linear. Our major result is that the tempera- 
ture Tp£) deduced from the (linear) FDR does not coin- 
cide with the statistical temperature Tth, and that Tpo 
does not take equal values in two subsystems when one 
considers an inhomogenous version of the model. This 
suggests that FDR are not necessarily the relevant way 
to define a temperature in the context of non glassy out- 
of-equilibrium steady-state systems. 

In addition, a numerical renormalization procedure 
suggests that detailed balance is generically restored on a 
coarse-grained level when it is not satisfied by the micro- 
scopic dynamics. This renormalization procedure yields 
a new parameter rj,, describing the deviation from equi- 
librium, which can be analytically computed within a 
mean-field approximation. This leads to a macroscopic 
description of the system with two parameters, namely 
Tth and 1]^. 

Finally, from a more general point of view, the present 
work raises important questions concerning the way to 
extend the concepts of statistical mechanics and ther- 
modynamics to out-of-equilibrium systems. On the 
one hand, the very definition of thermometers in non- 
equilibrium systems appears to be a highly non triv- 
ial issue, as the way to couple the thermometer to the 
system is not unique. Thus one may need to impose 
some -still unknown- prescriptions on the coupling to 
get a well-defined measurement. On the other hand, the 
present work may be of some relevance for the descrip- 
tion of non-equilibrium systems in which a global quan- 
tity is conserved. For instance, one may think of the 
two-dimensional turbulence where the vorticity is glob- 
ally conserved [s^, 0, 113 , or of dense granular matter 
in a container with fixed volume, in which the sum of 
the local free volumes would also be conserved. Indeed, 
the present model, for which the probability distribution 
is generically non uniform over the mutually accessible 
states (i.e. states with the same value of the energy, or 
volume, etc.) may allow in particular to go beyond the 
so-called Edwards' hypotheses 0,11^113, according to 
which all accessible blocked states have the same proba- 
bility to be occupied. 
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APPENDIX A: TIME-DEPENDENT ENTROPY 

In this appendix, we show that the time-dependent en- 
tropy Ssit) defined in Eq. is a non-decreasing func- 
tion of time. Let us first recall its definition: 



N 

J Y[dx,P{{x,},t) In 



i=l 
iN 



PEi{x^},t) 



(Al) 



with /({x,}) = nr=i \g'{x^)\gix,)'^-\ and PE{{x,},t) is 
the probability measure restricted to the hypersurface of 
given energy E: 

P{{x,},t) = PEi{x^},t) S (^^9ix^)-E^ (A2) 



Taking the derivative with respect to time, one finds: 

N 

(A3) 



cLSe 



dt 



dP , 
'dt 



Yl^dxi ■^i{xi},t) In 



PE{{x^},t) 



f{{x.}) 



since the integral of the time derivative of the logarithm 
vanishes. One can then use the master equation to ex- 
press dP/dt as a function of P{{xi}, t) and of the transi- 
tion rates. The obtained expression can be symmetrized 
by permuting the integration variables Xi and to get: 

^ = 2 /n^^'^^^ [W{{x'^\{x,})P{{x.},t)~ 

i—1 

Wi{xMx'.})P{{x'^,t)] X 

, PE{{x^},t) PEi{x',},t) ] 

f{{xi\) .f{{x[\) \ 

Now one can use the detailed balance relation Eq. H12|l : 
W{{x[]\{x,]) /({xj) = W{{x,}\{x[]) f{{x[]) (A5) 
and write dS/dt in the following way: 

. N 

l[dx^dx',W{{x'Mx,}) f{{x^}) X 



dS 
'dt 



i=l 



PEi{x^},t) PEi{x'i},t) 

fiiXr}) f{{x'J) 
^^ PE{{x.},t) ^^^PE{{x'A,t) 



(A6) 



fi{x^}) /(Kl) 

In this form, it is clear that the time derivative of the 
entropy is always positive. It vanishes only for the steady 
state distribution: 

PE{{X^}) 



1 



Zn{E) 



.f{{x^}) 



(A7) 



and the corresponding maximum value of the entropy is 
equal to: 



S{E) = \nZN{E) 



(A8) 
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APPENDIX B: MEAN-FIELD MASTER 
EQUATION 



In section HV CI a simple steady-state master equation 
was introduced to describe the one-site distribution p{£i) 
of the energy Si = g{xi) -see Eq. in the case of 

an arbitrary distribution ili{q). We show here how this 
simple equation can be derived from the master equation 
associated to a A^-site model with infinite range interac- 
tions. Introducing such long range interactions is a usual 
way to build a mean-field version of a model. To be more 
specific, we generalize the model introduced in Eq. Q in 
order to allow redistributions over any pair of sites (j, fc), 
and not only on the links of the lattice. As a result, the 
lattice becomes useless in this version of the model. 

The transition rates read: 



j<k 



(Bl) 



6{e'j + e'fe - Ej - £k) / dqij{q) 5 [e'^ - qisj + e^)) 
Jo 

where the sum runs over all pairs (j, k) . The factor 1/A^ is 
introduced so that each site keeps, in the thermodynamic 
limit ^ oo, a probability per unit time of the order of 
one to be involved in a redistribution. 

The stationary distribution PMF{{£i}) satisfies the fol- 
lowing master equation: 

N 

PmfHs.}) l[d<W{{e'Ms.})= (B2) 

•' i=l 

. N 

l[de',W{{e..}m)PMFm) 

i=l 

The first integral is the total exit rate from configuration 
{ej, and is equal to (TV - l)/2 from Eq. ||BT|) . So the 
last equation can be rewritten in a more explicit form: 



PMF{{£^}) 



N{N - 1) 



E 

j<k 

r-l 



£k) X (B3) 



dqip{q) 6 {ej - q{E'^ + 4)) Pmf(£j-, 4: {e«}«#j,fc) 



In order to go further, one has to assume that the distri- 
bution PMF{{£i}) factorizes: 



N 



PMFi{e^})=l[piE,) 



(B4) 



where p(e) is the one-site distribution. This assumption 
is justified in the limit of large N. Integrating over all 



variables except ei, one gets: 



JV-l N 



^ N(N-l) 
=1 fc=,+i ^ ^ 



(B5) 



Y[dE^ Yl p{E^) dE'jP{E'j) dE'kP{E'k)x 

»=2 ^=iJ,k •'0 Jo 

^i^'j + e'fe - £j - £k) [ dq ilj{q) 5 {ej ~ q{E'j + e'fe)) 



The r.h.s. can then be decomposed into two terms, one 
corresponding to z = 1 and the other one to i > 1, which 
are called respectively Ri and R2 in the following: 



p{ei) =Ri+R2 



(B6) 



The first term Ri is associated with redistributions in- 
volving site J = 1 as well as another arbitrary site k. It 
is actually independent of k, so that i?i is the sum of 
{N — 1) identical terms. Integrating over Ek removes the 
delta distribution S{e'i + e'i. — Ei — Ek), and one finds: 

R, ^ - dE[p{E[) dE'^piE'^) X (B7) 



dq'4}{q)5{Ei - ^(e'l +£2)) 



On the other hand, the second term R2 is the contri- 
bution from all the redistributions involving sites j = 
2 ... TV, but not site j = 1. There are {N - 1){N - 2)/2 
such pairs of links, which all give the same contribution 
to i?2- So i?2 can be written: 



N — 

i?2 = -^P{£l) I dE'2p{E'2) I dE',p{E',) X (B8) 



dq / dE2 S {e2 - qiE'2 + £3)) 



N 



All the integrals in the r.h.s. of the above equation give 
a contribution equal to unity, so that R2 reduces to: 



i?2 = ( 1 - ^ ) P{ei) 



(B9) 



Replacing the above results into Eq. (|B6|I . one finally 
obtains the following equation: 



Pie) 



nOC nOC 

/ deip(ei) / dE2p{E2) X (BIO) 
Jo Jo 



dqij}{q) (5 (e - q{Ei + E2)) 



which is precisely Eq. (|85|l . 
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with equal probability. 

Actually, only the behavior of g{x) in the vicinity of a; = 
is responsible for the divergence of the susceptibility x(*)- 
For an even regular function g{x) such that g"{0) 7^ 0, 
one has g{x) ~ for a; ^ 0, and the response remains 
linear in h for all positive value of ry. 
The same conclusions hold for more general functions 
g{x), but the results then take a less concise form. 



